RU06-06-B 



Low Shear Viscosity due to Anderson Localization 
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We study the Anderson Localization effect on the shear viscosity in a system with random medium 
by Kubo formula. We show that this effect can reduce the shear viscosity nonperturbatively. Then 
we discuss its possible implementation in heavy-ion collisions, where the created heavy bound states 
or other collective modes may play the role of the random scatterer underlying Anderson Localization 
effect. 
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One of intriguing properties regarding the newly iden- 
tified quark-gluon plasma (QGP) created by relativistic 
heavy ion collisions [J, 0, H, 0] is its near perfect hydro- 
dynamical behavior. In particular the ratio of the shear 
viscosity to the entropy density, deduced from the ellip- 
tic flow measurement is very small ry/s ~ 0.2 [5]. The 
weak coupling calculation gives a larger ratio of order 
5.12(<7 4 In for g « 1 for three flavors of massless 
quark 0, |r| > with g the QCD running coupling con- 
stant, which becomes unreliable near T c where g > 1. 
Non-perturbative effects are expected to play an instru- 
mental role to explain the observed ratio. 

An analytical calculation of QCD beyond perturba- 
tion theory is not available and the lattice simulation of 
transport coefficients [l8[ suffers from large error bars as- 
sociated to the analytic continuation to the mass shell, 
and so far only for pure SU(3) Yang-Mills theory with- 
out quarks. A lower bound of the viscosity-entropy ratio 
was estimated using the uncertainty principle within the 
framework of kinetic theories ^ > y^Q- This bound 
is surprisingly close to that of the J\f = 4 supersym- 
metric Yang-Mills theory at large N c and large 't Hooft 
coupling following the conjectured AdS/CFT duality Q 
2. = i ~ 0.08. The experimental result is rather close 
to these bounds and the underlying QGP seems strongly 
interacting 

Besides the lower bound of the viscosity-entropy ratio, 
there is little understanding of the physical mechanism 
that contributes to the suppression of the shear viscos- 
ity. One possibility is the existence of many zero energy 
bound states as was suggested in Q . Another possibility 
is the anomalous viscosity in a weakly but an expanding 
QGP within the formalism of Boltzmann equationfioj]. 

In this letter, we shall explore a new mechanism that 
may lead to a small viscosity in a system dominated 
by elastic scattering. This is the Anderson Localization 
(AL) effect studied in the context of condensed matter 
physics. 

The prototype AL effect refers to the multiple scatter- 
ing process of a wave propagating in a disordered medium 
where the individual scattering is elastic. Naively, one ex- 
pects that the net scattered wave equals to the incoherent 
sum of individual ones but a closer look reals that this 



is not the case. The scattering amplitude following any 
path of multiple scattering is in phase with the one fol- 
lowing its time reversed path when the scattering angle 
reaches 180, independent of the location of the individ- 
ual scattering events along the path. The elasticity or 
approximate one is required to maintain the phase rela- 
tion between the two scattering waves. This coherence 



effect was first suggested by Anderson [12j and was veri- 
fied experimentally by the intensity peak of the reflected 
light off the surface of an amorphous material with nor- 
mal incidence [l3} (the figure 1 of this paper provided an 
intuitive illustration of the AL effect). The field theo- 
retic treatments of AL have been developed in the con- 
text of the electrical conductivity [3] in an amorphous 
metal and the energy transport of light wave in a random 
medium 

In what follows, we shall generalize the existing analy- 
sis for the Anderson Localization effect on the electrical 
conductivity and the energy diffusivity to the momen- 
tum transport coefficient, shear viscosity, by considering 
a massless scalar field in a random medium. The mass- 
less particles resembles the thermal partons in the QGP 
produced by RHIC. Then we assume that there are some 
random scatterers required for the AL effect, for instance, 
the heavy bound states produced in the heavy ion colli- 
sions, or reminiscent of collective modes of the initial 
state. Because of the crudeness of our model, our con- 
clusion remains qualitative. 

While the transport coefficients can be calculated by 
diagrammatic expansion of the exact Kubo formula, 
they are usually extracted from kinetic theory by the 
Chapman-Enskog approximation of the Boltzmann equa- 
tion. It has been shown that the latter approach is equiv- 
alent to the ladder re-summation of the former one for self 
interacting scalar field and pure Yang-Mills field 17, 
The coherence effect underlying the Anderson Localiza- 
tion corresponds to a set of maximally crossed diagrams 
which is beyond the ladder ones, and therefore can not 
be obtained from the kinetic theory. 

Let us consider a simple model of massless scalar field 
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in a random medium. The Lagrangian reads 
1 dd> d 



where 



C 



2dx fl dx^ 



K(f> 2 , 
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FIG. 1: . Diagrammatic representations 



Here u(r — Rj) is the potential of an impurity at the 
random position Rj . The physics of AL is not sensitive 
to the details of the coupling. The Kubo formula of the 
shear viscosity is 



v = 77T lim n l im lm x(q, w) 

II) uJ—>0 q— >0 



(3) 



where x{Qi w ) is the Fourier transformation of the re- 
tarded two point Green function of the traceless stress 



tensor, -Ki 



d4> d<j> 
dxi dx-j 



d 4 Xe-* Q X < Trp[7r y (X),^(0)] > 8(t). 

(4) 

where p — Ze & is the the density operator with 
Trp = 1 . The bracket < ... > stands for the aver- 
age over the centers of individual impurity potentials, 
Rj . The four vector notation X — (r , it) for coordinate 
and Q = (q, ito) for momentum-energy are adapted. The 
two point function in eq.Jl]) can be manipulated with 
the CTP( closed time path ) formulation [20( or by more 
traditional means [14|. We find 



X(<f,w) = 



d 3 p f d 3 p' 



hj{p,q)hj{0,-q) 



(2tt) 3 J (2tt) 3 ~ 
I" ^:{[n( P o)~n(p^ RA (P,P',Q) 
n (po ) $ RR (P, P', Q) + n (p+) t> AA (P, P', Q)} 



X HA (q, to) + X HH (q, w) + X AA (q, to) 



(5) 



where n(p ) = (e Ppo - 1) , p* = p a ± \u) and I lj (p, q) 
is the Fourier component of the derivative operators in 
7T;, . We have 



h 3 {p, 0)lij(g, 0) = -p 2 p' 2 P 2 (p ■ p') 



(0) 



with Piip-p') the second Legendre polynomial. The func- 
tion $ a p{P, P', Q) with a and (3 equal to R(retarded) or 
A(advanced) can be expressed in terms the retarded and 
advanced full boson propagators 
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R{A) 



(P) = 



pI -p 2 



(7) 



together with the 1PI vertex functions T a g(P, P' ,Q), i.e. 
$ a0 (P,P',Q) = G a (P+)G (P-) x (8) 
~{2*) 3 5 3 (p- f) - iT a0 (P,P' ,Q)G a {P' + )G {P'_) 
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FIG. 2: . Schwinger Dyson equation and it's diagrammatic 
expansion 



where T, R (a)(P) is the retarded(advanced) self-energy 
function, P± = P±\Q, P' = (f, ip ) and P" = (p" , ip ). 
We have = pt^(p ) + 0(p 2 ) where tp^pa) de- 

notes the T-matrix of the scattering by the potential 
u(f) and p is the density of the impurity centers. The 
Born approximation of T matrix reads tp'.p(po) = up>-p 
with upi-p the Fourier transformation of u(r). We have 
ImSfl(P) = —ppo&(po) + 0(p 2 ) with cr(p ) the cross- 
section of the scattering by a single term of 

Since only $>ra(P,P' ,Q) contributes to the pinching 
singularity of the po-hitegration at weak coupling, we 
shall focus on it and see how the AL effect suppresses 
its contribution. The subscript R, A of T will be sup- 
pressed. The diagrammatic representation of §ra to- 
gether with the elements of the diagrams are displayed 
in Fig.l, where the upper(lower) solid line stands for the 
retarded(advanced) 0-propagator and the dashed line de- 
notes the impurity induced interaction. The bare vertex 
of Fig.lc reads V = -2inptp> + tf + (p£)t%,_ f _(po)6(p' - 

po), from which the full vertex T is generated. The en- 
ergy delta function stems from the elasticity of the under- 
lying scattering, which leaves the energy running along 
each solid line of Qra conserved and all loop integrals 
inside $ra over spatial momenta only. It follows from 
the Dyson-Schwinger equation for r(P, P',Q), shown in 
Fig. 2, that 



T(P,P',Q) = f(P,P',Q)- 

j^T(P, P", Q)G R (P!l)G A (P'l)T(P", P', QJ9) 
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which reduces to 

2( Po u + 2i 1Pa -p- q)f(P, P', Q) = f (P, P', Q) 

where f(P, P' , Q) = T(P, P', Q)/[2(p u + 2i 1Po - p- q)}, 
AG(P,Q) = G R {P+) - G A (P~) and f(P,P',Q) is the 
2PI part of T(P, P',Q). The leading order of f reads 



r = -iptg f + (pt)tg_tf_{PQ )• 



(ii) 



For |£„(P) << p, AG(P, Q) is sharply peaked at p = po 
and we may approximate £t?(P) ~ — Ea(P) — — 2ipo7- 
Introducing the partial wave expansion 

T(PP\Q) = Y,(21 + ±MPiP',u)Pi(p-p') (12) 



at Q = (0,oj) and that of f(P,P',Q), an approximate 
solution to (HHJ) at small oj and q can be obtained 

r(P, P , Q) = ; . n 2 + 4tt7 2^ 7 P(p • P 



where 



li = 7 + 



uJ + iDq 2 ^ 



' / (^ f (P,P',Q)P(p.p')AG(P',Q) 



(13) 



4^o 7 (2?r) 



(14) 

with the diffusion constant D(ui) = and the mass 

shell approximation AG(P, Q) ~ 27r<5(p ! — p 2 )sign(po) has 
been employed. The emergence of the hydrodynamical 
pole ( or the absence of 70 ) is the consequence of the 
Ward identity Q 



X R (P)-E A (P) = 



(2tt) 3 



AG(P 1 ,Q)f(P 1 ,P,0). (15) 



Inserting (fTB")) into (8) and setting Q = (0,w), we find 

^ = ir 2 AG(P,Q)5(p-p>) ^ cq 
2jpq Attlo 

7 Ty (2l + l)Pi(p-pf) 



8in 



i>x 



7i 



The shear viscosity picks up I — 2 partial wave of 
$(P,P',Q). It follows from ©, ©, ©, OH) and the 
mass shell approximation that 



i(3 f°° 



^72(po)(e*o-l) 2 



(17) 



The sum of the ladder diagrams shown in the 2nd line 
of Fig. 2, that corresponds to the kinetic theory ap- 
proach, can be obtained from the general solution flB")) 
with T of (|12p and (fT4"]) replaced by its leading order (fTTj) . 



The corresponding 7; 's will be denoted by . We have 
cq = —82*777 and 



ladder 



0P,P',Q) = 



16l7T7 2 

— Doq 2 



a (°) 



(18) 



The result 



with the bare diffusion constant Dq 
of the shear viscosity following from the kinetic theory 
amounts to approximate 72 in the formula (|17p by 7^ ■ 
The AL effect induced by the backward coherent scat- 
tering is reflected in the maximally crossed diagrams 
shown in the third line of Fig. 2 which is a subset of 
2PI diagrams for T. The sum of this set of diagrams, 
U(P, P', Q) can be obtained from the that of ladder dia- 
grams shown in the first line of Fig. 2 by reversing one 
of the rails ( dropping the bare vertex ) . We have 



t/(p,p / ,Q)=r ladder (p reV) p r ' ev ,i 



(19) 



where P rev = 
and Q rev = (p 



P> 

rev 



\{p-p*+q),ipo 
+ pVw)03- We find that 

16i7T7 2 



1 1^ 
2 



U(P,P',Q) = 



D (p+j?) 2 



(j?-p+q),ipo 



(20) 



and f (P, P', Q) = f + U(P, P', Q) + .... Therefore the 
maximally crossed diagrams renders the integrand of (14) 
to diverge in the backward direction, p 1 — > — p, as u) — > 0, 
and this singularity extends to all partial waves. 

In space dimensions d < 3, the backward singularity is 
sufficient to make ji divergent at to = and thereby to 
make the dressed diffusion constant vanish there for ar- 
bitrarily weak scattering. But a critical strength of scat- 
tering is required for the purpose in d = 3. A self consis- 
tent treatment of the backward singularity by summing 
up the most singular set of maximally crossed diagrams 
amounts to replaced the bare diffusion constant in the ex- 
pression of U(P, P' ,Q) with the dressed one, D(u>) 21|. 
Eq. (fl4|) for I = 1 becomes then a self-consistent equation 
for D(u>). Note that both Dq and D(lo) depend on the en- 
ergy pq. Upon introducing a cutoff k c that restricts the 
momentum integration within the neighborhood of the 
backward scattering and approximating AG(P, Q) by its 
(16) peak value [221] . — , the self consistent equation for D(ui) 



becomes 
1 

D 



1 



24tt7 2 
Po 



d 3 f 



P-P 



(21) 

which is of the same form as that in citezbsu. In terms 
of the localization length defined in T5I], the limit 




(22) 



the onset of AL is characterized by £ ^ 0, which implies 
that lim^—yfj D(uS) — 0, and turns the hydrodynamical 
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modes into localized ones via iu — Dt^p + p 1 ) — > iu[l + 
(P '+ f) 2 /£ 2 ]- Carrying out the integration of (j2~Tj) . we 
obtain the equation for the localization length, 



fc c £ — tan 1 (k c (.) = 



12 7 ' 



A real and nonzero solution for £ exists for 



12 7 2 



< 1, 



(23) 



(24) 



where we have set the cutoff k c ~ 7. It follows from 
ea.([I3)l for I = 2 ( with D Q replaced by D in T ) that 
AL effect renders 72 (po) — > 00 when the condition (12"4"1> is 
met. While it is unlikely that the condition (f2"4")) holds for 
all po ( because of the unitary bound that prevents 7(^0) 
from growing indefinitely with po ) , a finite domain of po 
where AL takes place is quite possible. If the AL domain 
covers the peak of the function p^e^ / '{eP Vo — l) 2 in the 
integrad of the formula (fTT|) , a significant reduction of 
the shear viscosity is expected. In the parallel case of the 
Mie scattering of an electromagnetic wave by a random 
ensemble of metallic spheres studied in [l5[ , AL does oc- 
cur for po greater than a critical value. But the upper 
bound of the AL domain was not reported there. But 
we have to admit that the mass-shell approximation em- 
ployed to obtain the approximate solution (| 16[) becomes 
marginal for d = 3 because of the condition (1241) . 

The self-consistent equation parallel to (|2"Tj) for d = 2 
reveals that the localization length survives for an arbi- 
trarily week coupling and for all p Q . No further conditions 
like are required. 

In conclusion, we have proposed a non-perturbative 
mechanism, the Anderson Localization effect, that may 
suppresses considerably the shear viscosity and other 
transport coefficients. The AL effect is caused by multi- 
ple coherent scattering in the system, which goes beyond 
the Boltzmann equation. It can be treated systematically 
in d = 2 and self-consistently in d = 3. Although the AL 
effect is discussed in a simple toy model in this letter, 
the physics of AL is more general. It only requires the 
individual scattering event to be sufficiently elastic. Re- 
garding the QGP produced in RHIC, a potential source 
for AL is the scattering by the heavy bound states pro- 
duced in collisions, whose masses could be as heavy as 
2GeV. If such states behaves like NR particles, the mo- 
mentum transfer in a scattering event could exceed by far 
the energy transfer. Other non-equilibrium effect, such 
as the color gauge field produced in an expanding QGP 
[l(]| may also implement the elastic scattering necessary 
for AL. Much study is still needed to ascertain the rele- 
vance of AL to the low viscosity of sQGP and other RHIC 
phenomenology, say the energy losses process. 
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